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1 Introduction 



In the past few decades, the supersymmetric approach has been profitably apphed to 
many non-relativistic quantum mechanical systems [1-5]. The SUSYQM has provided sat- 
isfactory results concerning different non-relativistic quantum mechanical systems, such as 
the exactly solvable and partially solvable potentials [ 6-9] . The exactly solvable potentials 
can be understood in terms of a few basic ideas which include supersymmetric partner 
potentials, shape invariance and operator transformations. 

Among the interesting problems of the non-relativistic quantum mechanics which aim to 
find exact solutions to the Schrdinger equation for certain potentials of the physical interest, 
the screened coulomb potentials have been studied in a variety branches of physics such 
as atomic, nuclear and plasma physics [10-17]. Various types of the screened Coulomb po- 
tentials hke the Yukawa, Debye-Hiickel and exponential-cosine screened Coulomb (ECSC) 
potentials are discussed in non-relativistic quantum mechanics [18-20]. While the screened 
Coulomb potential which is in the vector coupling prescription leads to an exactly solvable 
for one dimensional Dirac and Klein-Gordon equations [21-23], the Schrodinger equation 
for these potentials is not exactly solvable. Perturbative and approximation methods have 
been applied to obtain their energy eigenvalues by using hypervirial/shifted 1/N expansion 
technique, variational approach, Pade approximates, numerical integration and group theo- 
retical approach [24-32]. The bound state energies of some potentials like the Morse, Poschl- 
Teller and other exponential type potentials are evaluated through the SUSYQM method 
by following the PT-symmetric formalism [30-36]. PT-symmetric Hamiltonians satisfy the 
parity (P) and time reversal (T). The eigenvalue spectra of PT-symmetric potentials may 
be real or complex [37]. If PP-symmetry is not spontaneously broken, the form of spectra 
is real. For a class of non-Hermitian Hamiltonians, the concept of pseudo-hermiticity is 
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valid [38]. In this work, the energy eigenvalues and corresponding eigenfunctions of PT- 
/non-PT-syminetric and non-Hermitian types of the Exponential- Cosine Screened Coulomb 
potential are obtained by using the Hamiltonian Hierarchy method within the context of 
PT-symmetric quantum mechanics (PT-SQM). 

This paper is organized as follows: In section 2, we give a brief pedagogical review of the 
Hamiltonian hierarchy method. In section 3, we apply this method for the Exponential- 
cosine screened Coulomb potential. In sections 4 and 5, the method is applied for the 
Pr-/non- PT-symmetric and non-Hermitian cases of this potential. In section 6, the re- 
sults are discussed as a conclusion. 



2 Hamiltonian hierarchy method 

The radial Schrodinger equation for some specific potential energies can only be solved 
analytically for the states with zero angular momentum [36, 37]. However, in supersymmet- 
ric quantum mechanics one can deal with the hierarchy problem by using effective potentials 
for non-zero angular momentum states in order to solve the Schrodinger equation analyt- 
ically. Hamiltonian hierarchy method suggests a hierarchy problem in the frame of the 
SUSYQM in which the adjacent members are the supersymmetric partners that share the 
same eigenvalue spectrum except for the missing ground state. 

In this method, the first step is to look for an effective potential similar to the origi- 
nal specific potential and inspired by the SUSYQM to propose a superpotential, namely 
VF(z+i)(a;), as an ansatz, where (/ + 1) denotes the partner number with Z = 0, 1, 2.... Sub- 
stituting the proposed superpotential into the Riccati equation, 

- El,,, = Wl,,,{x) - ^^llmM, (1) 
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the (/ + l)th member of the Hamiltonian hierarchy can be obtained. As a result, considering 
the shape invariance requirement [14] , the bound-state energies can be derived out through 
the Eq. (1), and the corresponding eigenfunctions by means of, 



(2) 



3 Exponential-cosine screened Coulomb potential 



The cosine screened Coulomb potential is written as. 



Substituting cos(/xr) 



V{r) = -^e-^"cos(Atr). 

in the above potential, we get. 



(3) 



(4) 



or, 



V(r) = -I 



{ifi-X)r _j_ g-(i//+A)r 



To simplify the calculations and for simplicity, let us take q = 2. Therefore, 



(5) 



„-(A-i/i)r- p-{\+iiJ,)r 

V{r) = -- . 

r r 

This potential can be considered as two separate parts as, 

p-{\-in)r 

V,{r) = , 



(6) 



(7) 



and. 



V2{r) = 



-{X+ilJ.)r 



(8) 
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By defining X — ijji = a and X + i/j, = P, the superpotential proposed as an ansatz for tlie 
Vi{r) potential becomes, 



(TP '^^ 1 (T 

w'W'-) = -(' + i)y3^ + 7^-2- P) 

According to tlie Hamiltonian liierarcliy metliod, tlie corresponding eigenfunction for tliis 
superpotential will be, 

(r) = (l - e-"'^)'^' e-(^-^)^ (10) 

Assuming that the radial trial wave function is given by (10), we replace a by the variational 
parameter /ii, and as a result, 

^^^(r) = (l - e-'^i'-)'^' e-^TT-i-'^y. (11) 
The variational energy is given by. 



E^ii — 



(12) 



The superpotential proposed as an ansatz for the V2(r) potential is, 

8e~^'' 1 3 
W,,^,,(r) = -(/ + !) + 7^ - f ' (13) 

and the corresponding eigenfunction for this superpotential becomes, 

^02(r) = (l - e-^')'^' e-^Th-^y. (14) 

Again, assuming that the radial trial wave function is given by (14), we can replace f3 by 
the variational parameter 112, and get. 
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^^^(r) = (l - e-'^^'-)'^' e-^Tk-'^y. (15) 
The variational energy is given by, 



(16) 



Thus, by minimizing the energies E^j^-^ and E^^ with respect to the variational parameter /^i 
and //2, one obtains the best estimate for the energy of the exponential-screened Coulomb 
potential. 

As the exponential-cosine screened Coulomb potential is not exactly solvable, the superpo- 
tentials given by Eqs.(9) and (13) do not satisfy the Riccati equation, but they do satisfy 
for effective potentials instead. Vie// and V2e// as. 



Vuff{r)^ ' ^ ' +E{fi^). (17) 

and 



V2eff{r)^ ' , ' +E{fi,), (18) 



where Wi — Wi{a — jii) and W2 — W2{P — 1^2) ■ A*i and 112 are the parameters that 
minimize the energy expectation values (12) and (16). They are given by, 

V,,„ (r) = + + ^' , + i( 1 - + E(a), (19) 

leffKj 1-e-"'- 2 (l-e-°'')2 2^/ + l 2^ ^ ^' ^ ^ 

and 



^^e// (r) = -3;^^ + 2 (1 - e-/^0^ 2^H^ " 2 ^ + ^^^^ 
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By substituting the values of a and /3 in (19) and (20), one can obtain the bound state 
energies of the exponential-cosine screened Coulomb potential as, 



1 A2-/x2 
+ 7-^ 



A 



(^ + 1) 



l + l 



(21) 



It is interesting to notice that for = 0, the exponential-cosine screened Coulomb potential 
reduces to the to the form called Yukawa potential and that Eq.(21) reduces to. 



.(^ + 1) 



A" 
2 



(22) 



4 Non-PT-symmetric and non-Hermit ian exponential- 
cosine screened Coulomb case 

The non-PT and non-Hermitian cosine screened Coulomb potential can be defined as. 



V{r)^ ^e-^''cos{iir). (23) 
r 

or simply, 

y(r) = |(^ + ^). (24) 
In this case the proposed superpotentials can be, 

and 
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iBe'P'^ 1 B 
W,,,,,{r) = + + - f ' (26) 

Though the superpotentials are complex, following the same method will yield the same 
energy eigenvalues as in (21). 



5 PT-symmetric and non-Hermit ian exponential-cosine 
screened Coulomb case 

The PT symmetric and non-Hermitian cosine screened Coulomb potential can be intro- 
duced as, 

T/(r) = -^e-^^"cos(//r). (27) 
r 



or, 



V{r) = 



^-i{X-lx)r _j_ g-i(A+//)r 



Taking, {X — ^) — — ckq and A + = /3o, we will have. 



„ p-iaor p-i0or 



and as a result the superpotentials can be proposed as. 



(28) 



(29) 



iyi(,+i)(r) = -(/ + !) 



1 - e-*"o'' I + 1 



and 



dne'^^°^ 1 



2 ' 



/3o 



1 - e-'*'- Z + 1 2 ' 



(30) 



(31) 



In conclusion, by applying the method the same energy eigenvalues will be obtained as in 
(21). 
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6 Conclusions and remarks 

We have applied the the Hamiltonian hierarchy method within the framework of the 
SUSYQM formulation by presenting a superpotential that yields a trial function to calculate 
the approximate bound state energies and corresponding eigenfunctions for the exponential- 
cosine screened Coulomb potential. We have also considered its different symmetric forms 
in our calculations. As the energy spectrum of the PT-invariant complex-valued non- 
Hermitian potentials may be real or complex depending on the parameters, we have clarified 
that there are some restrictions on the potential parameters for the bound states in PT- 
symmetric, or more generally, in non-Hermitian quantum mechanics. Furthermore, it is 
shown that the superpotentials, their superpartners and the corresponding ground state 
eigenfunctions satisfy the PT-symmetry condition. 

Finally, we can add that our approximate yet accurate results of complexified exponential- 
cosine screened Coulomb potential by the justification of the numerical results presented in 
Table 1 motivate an appropriate approach to analyze the exactly and non-exactly solvable 
potentials. We believe that this method may increase the number of applications in the 
study of different quantum systems. 
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Table 1 : Energy eigenvalues of ECSC as a function of the screening parameter A for I5, 
2p, 3rf and 4/ states in Rydberg units of energy. 



State Is 


Screening 


SUSYQM 


Hypervirial 


NR-QM 


Exact 


A 


Our Work 


Solution [40] 


Variational [20] 


Numerical [39] 


0.020 


-0.480290 


-0.480310 


-0.480300 


-0.480300 


0.050 


-0.451810 


-0.451800 


-0.451820 


-0.451800 


0.080 


-0.424560 


-0.424560 


-0.424570 




0.100 


-0.407070 


-0.407050 


-0.470600 


-0.407100 


State 2p 


0.020 


-0.211800 


-0.105890 


-0.211900 


-0.211900 


0.050 


-0.162500 


-0.080400 


-0.161500 




0.080 


-0.050500 


-0.046000 






0.100 


-0.092860 


-0.008000 


-0.092890 


-0.093070 


state 3d 


0.020 


-0.075020 


-0.037500 


-0.075030 


-0.075030 


0.050 


-0.033620 


-0.017340 


-0.033740 


-0.033830 


0.080 


-0.009020 


-0.008000 






0.100 


-0.038889 








state 4/ 


0.020 


-0.028750 


-0.014700 


-0.028970 




0.050 


-0.004100 


-0.003200 






0.080 


-0.184500 


-0.175000 






0.100 


-0.018700 
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